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In this paper, we introduce an asymptotic test procedure to assess 
the stability of volatilities and cross- volatilites of linear and nonlinear 
multivariate time series models. The test is very flexible as it can be 
applied, for example, to many of the multivariate GARCH models 
established in the literature, and also works well in the case of high 
dimensionality of the underlying data. Since it is nonparametric, the 
procedure avoids the difficulties associated with parametric model 
selection, model fitting and parameter estimation. We provide the 
theoretical foundation for the test and demonstrate its applicability 
via a simulation study and an analysis of financial data. Extensions 
to multiple changes and the case of infinite fourth moments are also 
discussed. 

1. Introduction. Univariate time series models, both linear and nonlin- 
ear, are as of today rather well understood. The literature furnishes a broad 
variety of contributions on the probabilistic structure of these models pro- 
viding, for example, criteria for the existence of stationary solutions and the 
finiteness of moments. Manifold estimation procedures are well-developed 
and tests assessing the structural stability and the goodness of fit are widely 
available (see, among many others, [9, 16, 18, 39, 40, 48] and the references 
therein) . 

This, however, is no longer the case for multivariate time series models, in 
particular for multivariate financial time series. While many generalizations 
to higher-dimensional settings have been proposed since the inception of the 
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univariate ARCH and GARCH models by Engle [23] and Bollerslev [11] (see, 
for example, the survey articles [7, 47] and Section 4 below), probabilistic 
and statistical tools are still in developing stages. The underlying nontrivial 
theory is often times only rudimentary developed. Establishing efficient es- 
timation procedures proves to be a challenging problem, since it bears the 
difficulty of incorporating a potentially large number of parameters into the 
algorithms. Tests that assess the structural stability of volatilities or cross- 
volatilities for multivariate nonlinear time series have yet to be introduced. 
On the other hand, multivariate nonlinear time series modeling is of high im- 
portance in practice as it is essential for an understanding of the coherence of 
exchange rates or market indexes. How movements of the financial markets 
are interpreted will significantly determine, for example, the asset allocation 
in a portfolio or the decision making of a risk manager (see [7, 24, 46, 47]). 
It is therefore crucial to understand the dynamic dependence structure in 
multivariate financial time series. 

The main aim of this paper is to add to the research in this area by 
studying in detail the volatility and co-volatility structure of <i-dimensional 
random vectors, which allow both linear and nonlinear specifications. We 
introduce an asymptotic test that is very flexible in that it only requires 
general but easily verifiable dependence assumptions to be satisfied by the 
underlying sequence. The procedure bypasses the difficulty of estimating 
a large number of model parameters and is very applicable also in high 
dimensions. As a byproduct, we provide stationarity conditions for several 
multivariate nonlinear time series models not previously established in the 
literature. 

The results are based on an approximation of the underlying random vec- 
tors with a sequence of m-dependent random vectors which in turn yields a 
multivariate functional central limit theorem. This approximation general- 
izes similar univariate results provided in [10, 32] and is the main theoretical 
tool in the statistical analysis. In practice, this approach is generally more 
convenient to apply than competing methods such as the various notions of 
mixing (see, e.g., [17, 43]) as it typically requires less restrictive and eas- 
ier verifiable assumptions. Moreover, the specifications employed here can 
even be modified to detect structural breaks in the co-movements of the d 
coordinate processes under the assumptions of heavy tails (less than four 
finite moments) and even an infinite variance-covariance structure. This is 
especially important in applications related to financial time series. 

An outline of the paper can be given as follows. In Section 2, we dis- 
cuss the model assumptions and state the main results. These results are 
then extended to multiple change scenarios and to sequences with infinite 
variance-covariance structure in Section 3. Section 4 deals with a number 
of examples that are included in the provided framework. Section 5 is con- 
cerned with practical aspects and investigates the finite sample behavior 
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through a simulation study and an application to stock index data. The 
multivariate functional central limit theorem is presented in Appendix A, 
while the mathematical proofs are given in Appendix B. 

2. Main results. Let (v.,- : j € Z) be a sequence of d-dimensional random 
vectors with i£[yj] = and .E[|yj| 2 ] < oo. Here | • | denotes the Euclidean 
norm in R d . The main aim of the present paper is to introduce a test- 
ing procedure that identifies structural breaks in the volatilities and cross- 
volatilities of the process (y,- :j G Z) based on observations of the random 
vectors yi, . . - ,y n . To do so, we consider the null hypothesis 

(2.1) tfo:Cov( yi ) = --- = Cov(y n ), 

which indicates the constancy of the covariances in the observation period. 
It is usually under Hq that a statistician can produce meaningful estimates 
and reliable forecasts. Whether or not Hq holds should therefore precede any 
further statistical analysis. It should, however, be noted that a common ap- 
proach to forecasting in the finance sector and elsewhere is to model volatili- 
ties via exponential smoothing with little regard to stationarity assumptions 
on the underlying processes. This works well as long as the volatilities can 
be assumed to evolve in a roughly smooth fashion, but leads to inaccurate 
forecasts in the nonsmooth break setting we consider in this paper. 

As the alternative hypothesis we specify in particular the common sce- 
nario which allows for at most one change in the covariances. More precisely, 
it is assumed that there exists a — typically unknown — time lag k*, referred 
to as the change-point, such that 

(2.2) H A : Cov( yi ) = • • • = Cov(y fc ,) + Cov(y fc , +1 ) = • • • = Cov(y n ). 

Generalizations to more general alternatives allowing several change-points 
are discussed briefly in Section 3 below. Since we are interested only in the 
second-order behavior of the random vectors (yj'-j 6Z), it is throughout 
this paper assumed that the expected values -E[vj] do not change over time. 
This may require an initial test of constancy of the means and, if necessary, 
a transformation of the data so that the expected values can be regarded as 
stable for the whole observation period. For further information on testing 
for changes in the mean, we refer to [18]. 

Even in the univariate case, there are relatively few contributions address- 
ing changes in the variance. Gombay, Horvath and Huskova [26] and Inclan 
and Tiao [33] discuss the detection of (multiple) changes in the variance of 
independent observations based on weighted test statistics and an iterated 
cumulative sum of squares algorithm, respectively. The special case of (inde- 
pendent) Gamma distributions has been treated in [31]. More recently, [37] 
have introduced tests for a single change in the volatility of an ARCH(oo) 
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series, while [20] analyze parametric nonlinear time series models by means 
of minimum description length procedures. There are, however, no in-depth 
treatments of (nonparametric) break detection methods for changes in the 
volatilities and cross-volatilities for multivariate, potentially nonlinear time 
series. 

To construct a test statistic for distinguishing between Hq and Ha, we 
let vech(-) be the operator that stacks the columns below the diagonal of a 
symmetric dx d matrix as a vector with (3 = d{d+ 1) /2 components. Assume 
for the moment that E[yj] = 0. Then, it holds under Hq that the expected 
values E[vech(yjyJ)] are the same for j = 1, . . . ,n. Consequently, a version 
of the traditional cumulative sum (CUSUM) statistic can be constructed 
using the quantities 



1 / ^ k n 

— £ vech[y,yj] - - £ vech^yj 

V \j = i j = i 



,n, 



which basically compare the estimators of E[vech(yjyJ)] based on k ob- 
servation with the one computed from all observations. For j = 1, . . . ,n, let 
Yj = Yj ~ Yn with y n = £ J2]=i Yj- If E[y ] ¥= 0> then the S k can be replaced 
with the mean corrected modifications 

1 ( k k n \ 
S k = -p ^ vechfy^yj] ^ vechfy^yj] , k = l,...,n. 

V n \j=l n j=l ) 

Before we introduce the precise form of the test statistic, we detail the 
assumptions needed on the sequence (y^- :j £ Z) . We are interested in several 
competing parametric specifications which are described in detail in Section 
4 below. It is, however, often times hard to distinguish between different 
parametric models in practice and therefore assumptions imposed on (yj : j £ 
Z) need to be general. Let a > 1. For a random vector x £ M. d , let |x| be its 
Euclidean norm in M d and ||x|| Q = (-Eflxl"]) 1 /" be its L a -norm. 

Assumption 2.1. Let (yj-j £ Z) be such that, under H , it satisfies 
the relations 

(2.3) y j = i(e j ,e j - 1 ,...), j £ Z, 

where f :M rf ' xo ° — > M d is a measurable function and (e,- :j £ Z) a sequence 
of independent, identically distributed random variables with values in M d> . 
It is further required that, under H , there is a sequence of m-dependent 
random vectors (y^ : J £ Z) such that 

(2.4) y ; m) =f( m )(£,,..., £j - m ), jez, 
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with measurable functions f( m ) : ]g>d'x(m+i) _> anc [ 

(2.5) Ellyo-yf ) m<oo. 

m>l 

Condition (2.3) states in other words that (yj-j £ Z) admits a causal 
representation, possibly nonlinear, in terms of the d'-dimensional innovation 
sequence (ej :j £ Z). A weak dependence structure is enabled through the 

introduction of the m-dependent random vectors (y^ in (2.4) which 

are close to the original (yj :j £ Z) in the sense of the distance measure (2.5). 
Note, however, that (y, : j £ Z) itself need not be m-dependent. This type of 
condition has first been used in the context of univariate random variables in 
[32] and later also in [10, 29]. We will show in Section 4 that several popular 
multivariate time series models satisfy Assumption 2.1. 

As shown in Appendix A, Assumption 2.1 induces a functional central 
limit theorem (FCLT) for the sequence (vechfy^yj] :j £ Z). FCLTs are fre- 
quently used in the asymptotic theory of weakly dependent processes. Var- 
ious forms based, for example, on near epoch dependence and mixing con- 
ditions are discussed in [19, 43, 45]. The advantage of Assumption 2.1 over 
its mixing competitors is twofold. First, it is tailor-made for univariate and 
multivariate GARCH-type processes as the construction of (y^ : j € Z) is 
virtually always straightforward (see [29] and Section 4 below), whereas 
mixing conditions are typically hard to verify and not applicable for several 
important time series models (see Section 4 below and [3, 29]). Second, mix- 
ing conditions require usually additional smoothness or restrictive moment 
assumptions (see [17]), while Assumption 2.1 gets by with finite fourth mo- 
ments (which cannot be improved on when studying the covariance struc- 
ture as in the present paper). These issues have led several authors such 
as [22, 50] to introduce more general dependence concepts which are more 
widely applicable. It is not immediately clear, however, whether the sharp 
almost sure invariance principles in [50] can be extended with reasonable 
effort to a multivariate setting. The methods in [22] are very general, but 
require more restrictive moment assumptions. In Section 3, we will address 
how to further generalize our methodology to multivariate sequences with 
infinite moments of less than fourth order. 

Condition (2.5) yields in particular (see Theorem A. 2 in Appendix A) 
that the long-run covariance 

T, = J2 Cov(vech[y yo ], vech[y jy J]) 

jez 

converges (coordinatewise) absolutely. Statistical inference requires the es- 
timation of S. We assume hence that there is an estimator S n satisfying 

(2.6) |E n -E| B = op(l) (n— ►oo), 
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where, for a d x d matrix M, \M\e = sup^g l-W^I/l^l denotes the matrix 
norm induced by the Euclidean norm on M. d . We shall discuss specific estima- 
tors S n obeying the rate in (2.6) in Section 5. We now suggest the following 
two sets of test statistics to test Ho against Ha- Define 

1 n 

A n = max S^t^Sk and fi n = - V SjfS" 1 ^, 

l<fc<n n — ' 

h=\ 

as well as 

1 n 

A n = max S^t^Sk and O n = - V SjS" 1 ^. 

Kk<n n 

k=l 

The limit distributions of all four test statistics are given in the next theorem. 



Theorem 2.1 (Asymptotic under Hq). Suppose that (y~j :j GZ) are d- 
dimensional random vectors satisfying Assumption 2.1 and let S ra be such 
that (2.6) holds. Given that E\yj] = 0, it holds under (2.1) that 



(2.7) A n A(a) = sup B e(t) ( n 



oo 



0<t<l 



1=1 



and 

(2.8) n n ^n{d) = J2 Bj(t)dt (rwoo), 

where d = d(d+ l)/2, (B^(t) :t E [0, 1]), 1 < I < D , are independent standard 

T> 

Brownian bridges and — > indicates convergence in distribution. The limit 
results in (2.7) and (2.8) remain valid if A n and fi n are used in place of A n 
and Q n , respectively. 

The proof of Theorem 2.1 is relegated to Appendix B.l. The distribution 
of the random variables A(0) and £l(d) were derived in [36]. It is shown there 
that, for x > and > 2, the distribution function of f2(0) has the series 
expansion 

p(nco) < x) 

(2.9) 

2 (0+1)/2 ^ r(j+V2) f , +0/4 , 2/ / 2j + 0/2 

j!T(D/2) °('-2)/2^ x i/ 2 

where V denotes the Gamma function and C( _2)/2 the parabolic cylinder 
functions (see, e.g., page 246 in [25]). Somewhat more involved formulas for 
the distribution of A(t>) are also provided in [36]. Since t) could potentially 
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be large, it is interesting to consider the distributions of A(5) and 17(0) as 
d — > oo. We therefore define the standardized variables 

j. . ACO) -Hi A , = , , -0/6 



Remark 2.1. 7/0 -> cx), i/ien both A(d) -^iV(0,l) andn(d) iV(0, 1), 
where N(0, 1) denotes a standard normal random variable. 

The proof of Remark 2.1 may be found in Appendix B.l. Next, we turn 
our attention to the behavior of the test statistics under the alternative 
hypothesis Ha- To establish a consistency result, we need another set of 
assumptions which will be specified now. 

Assumption 2.2. Let k* = [6n\ with some 6 G (0,1). Let {yj-j G Z) 
be such that the pre-change sequence (y,- :j < k*) is strictly stationary and 
ergodic with ^[yo] = and E^jyol 2 ] < oo, and the post-change sequence 
(y,- :j > k*) satisfies the relations 

(2.10) yi = y* + Zj> , j>k* + l, 

where (y* : j G Z) denotes a strictly stationary, ergodic sequence with E[yq] = 
and -Bflygl 2 ] < oo, and (zj jn : j G Z) is such that 

(2.11) max |z,-„| = op(l) (n— >oo). 

k*<j<n ' 

As common in the literature, the time of change k* depends on the sample 
size n (see [18]). This implies that we are actually dealing with an array 
under the alternative. For the sake of simplicity, this is suppressed in the 
notation. Roughly speaking, however, we are observing a stationary sequence 
until k*. After the change, a new sequence starts with yk*+i as initial value. 
By assumption (2.11), this sequence is asymptotically stationary. The same 
condition (2.11) also ensures that the law of large numbers applies to the 
post-change sequence. 

If the corresponding model parameter values change at time lag k*, the 
conditions imposed by Assumption 2.2 are, in particular, satisfied for all of 
the specific parametric time series models discussed in Section 4 below. The 
following theorem establishes the consistency of the test procedures. 

Theorem 2.2 (Asymptotic under Ha)- Suppose that (yj :j G Z) are d- 
dimensional random vectors satisfying Assumption 2.2. Then it holds under 
(2.2) that 

sup -^=S [nt ^ -S*{t) 
te[o,i] V n 



op(l) (rwoo), 
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where 

<7*m = /*(!- ^)^[vech[y y^ - vech[yS(y5) T ]], t G [0,0], 
W I «(1 - i)^[vech[y y^] - vechfy^f]], t G (0, 1]. 

suitable modifications of the limit process (S*(t):t G [0,1]), a corre- 
sponding statement holds true if the process (5[niJ '-t £ [0)1]) * s replaced with 
(5 LntJ :t€[0,l]). 

Observe that if the mean of the vech[y.,yj] changes, then the process 

{S\nt\l \fn'-t ^ [0,1]) will have a nonzero limit with probability one. This 
shows that our test has asymptotic power one. The proof of Theorem 2.2 
is relegated to Appendix B.l. Once Hq is rejected, a practitioner needs to 
locate the break point k* or, alternatively, the break point fraction (relative 
to the sample size) 0. For this, we suggest the estimator 

(2.12) n = - argmaxSjS" 1 ^. 

n l<k<n 

From Theorem 2.2, one can deduce intuitively that for large enough sample 
sizes n, the quadratic form <Sj^j ^^"^[nt] wm reach its maximum in the same 
argument as the limiting quadratic form (S*(t)) T T>~ l S*(t), thus yielding 
consistency. In fact, strong consistency of 9 n can be proved under quite 
general conditions. Although a detailed discussion of this is beyond the scope 
of this paper, we will devote several paragraphs in Appendix B.2 to show 
that the relation 

(2.13) \§ n -e\ = o( l ° gl ° gn ^ a.s. 

holds in many practical situations. 

We continue with extensions of our methodology in Section 3 and exam- 
ples in Section 4, while Section 5 examines practical aspects of the testing 
procedures via a simulation study and an application. 

3. Extensions. We discuss in this section two directions of possible ex- 
tensions of the theory introduced in Section 2. First, we show how our results 
may be applied also in the case of multiple changes in the underlying data. 
Second, we extend the approach to sequences of random vectors with heavy 
tails. 



3.1. Multiple break point detection. In the literature of break detection 
schemes for univariate processes, a number of multiple change-point testing 
and estimation procedures are well established. They are, however, mainly 
concerned with mean changes (such as [6]), or require parametric model 
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assumptions (such as [20]). The references for the univariate change in the 
variance case listed in the previous section have similar drawbacks. There 
are two papers by Andreou and Ghysels [1, 2] which address testing for 
multiple changes in the co-movements of financial assets from an empirical 
point of view without pursuing the corresponding asymptotic theory. This 
paper is thus novel in this direction. To incorporate the case of multiple 
breaks in the covariances into our setting, we can exchange H A with the 
following counterpart. Assume that, for some r > 1, there are time lags 
1 < k\ < ■ ■ ■ < k* < n such that, for s = 1, . . . , r, 

H* A : Cov(y fc *_ i+1 ) = • • • = Cov(y fc «) / Cov(y fe » +1 ) = • • • 

= Cov(y fc:+i ), 

where k$ = and k* +1 = n. If the change locations have the standard form 
used in the analysis of multivariate changes and if the change size is either 
constant or shrinking at an appropriate rate (see [6, 18]), then, all four test 
statistics A n , £l n , A n and O n have asymptotic power one. To make this state- 
ment precise, we replace Assumption 2.2 with the following requirements. 



Assumption 3.1. Let k* = [nd s \, s = 0, . .. ,r + 1, with = 6> < 1 < 
■ ■■ <9 r < 6 r+ \ = 1. Let the sequence (y^ : j = 1, . . . ,n) be such that y^ = y^ , 
j = 1, . . . , k\ and 

Yj — Yj,s ' z j,ni J — K s ? • • • ) K s+l ~ x > s — x ; ■ • ■ i r i 

where, for s = 0,...,r, (y* s :j£Z) denotes a strictly stationary, ergodic 

sequence with -E[yo jS ] = and -E[|yo, s | 2 ] < oo, and, for s = 1, . . . , r, {zf^ :j G 
Z) is such that 

max max \z^„\ = op(1) (n^-oo). 

l<s<rfc s *<j<fc s * +1 ' 

With Assumption 3.1, we obtain the following counterpart of Theorem 
2.2 under the more general alternative H A . A simple computation shows 
that the one-change case is reproduced if r = 1 is chosen. 



Theorem 3.1 (Asymptotic under H A ). Suppose that (yj :j = 1, . . . , n) 
are al- dimensional random vectors satisfying Assumption 3.1. Then it holds 
under H A that 



sup 

te[o,i] 



op(l) (n->oo), 
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where S*(0) = and 

r-l 

S* r (t) = £ " e u )E[vech[yl u (yl u ) T ]] + (t- # s )£[vech[y^(y^) T ]] - 1& 

u=0 

for t £ (9 S ,9 S+1 ], s = 0,...,r, with 

^ = ^(0 s+1 -^)E[vech[y o * s (y^f]]. 

With suitable modifications of the limit process (S*(t):t € [0,1]), a corre- 
sponding statement holds true if the process (S\ nt \ :t£ [0,1]) is replaced with 

(5 LntJ :t€[0,l]). 

One can now easily verify that sup te r 0) i] |«5*(<)| > under H\, so that all 
four test statistics have asymptotic power one. In the empirical Section 5, we 
illustrate that the binary segmentation procedure based on Vl n produces very 
reasonable results in an application to financial data and is therefore capable 
to detect multiple breaks in the variance-covariance structure of multivari- 
ate financial time series. Moreover, applying the technique developed in [6] 
would also imply the weak consistency of the estimators of 9i,...,0 r , which 
are defined analogously as in (2.12) and are evaluated utilizing the binary 
segmentation procedure. 

3.2. Transformations for sequences with less than four finite moments. 
There is plenty of evidence in the literature that the assumption of four 
finite moments for financial time series is rather restrictive in applications. 
Theoretical results such as limit theorems for the sample autocorrelation 
function and the quasi-likelihood estimation of the unknown model param- 
eters has for heavy-tailed univariate ARCH and GARCH processes been 
discussed in [21, 27, 41], among others. Even in the univariate case, these 
results are generally difficult to apply in practice because the statistical infer- 
ence procedures require, for example, precise knowledge of the observations' 
number of finite moments. We are unaware of any results concerned with 
the corresponding multivariate versions of heavy-tailed GARCH processes. 
The approach introduced in Section 2, on the other hand, is in principal not 
hampered by these limitations as the following argument shows. Let again 
(yj : i ^ ^) denote the multivariate process to be observed. For 5 £ (0, 1], let 
(xj :jeZ) be the fractional process transformation given by 

(3.1) x j = \y j \ 5 = (\y j , 1 \ 5 ,...,\y j>d \ s ) T , j G Z. 

Instead of detecting breaks in the covariance structure of the original (y,- : j £ 
Z), the transformed process (x-j'j £ Z) can be utilized. This bears the 
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further advantage of increased flexibility of the testing procedure as it is 
in the transformed version also possible to find, for example, breaks in 
the marginal distribution of the yj. To apply our theory to (3.1), condi- 
tion (2.5) needs to be verified for the Xj. Let therefore = |y^ I* 5 be 
the sequence of m-dependent variables used for the approximation. Then 
|| Xj — x^. m -' 1 1 4 < d^Hyj — yj- H45 and consequently 

(3.2) £ ||xo - IU < rf 1/4 E llyo - yt ] Wis < oo 

m>l m>l 

provides us with a sufficient criterion. Picking 5 < 1/2 shows moreover that 
our method works theoretically even if the process (y, : j G Z) does not pos- 
sess a finite variance-covariance structure. In the case of GARCH sequences, 
for example, it appears convenient to pick 5 = 1/2. In this case, (3.2) is even 
easier to verify than the original version in (2.5). We further discuss the 
choice 5 < 1/2 for the constant conditional correlation GARCH model in 
Section 4.2 below. It would be worthwhile developing a data-driven pro- 
cedure to determine an optimal value of 5. Heavy-tailed random vectors 
naturally display a greater amount of variability compared to light-tailed 
ones, so that it is generally much harder to discriminate between wild but 
stationary fluctuations and an actual break in the underlying process struc- 
ture. On the other hand, choosing 5 < 1 usually flattens the signal, so that 
breaks may be harder to detect. A selection procedure for 5 needs to balance 
these two properties. We leave this issue to future research. 

4. Examples. In this section, we discuss the leading examples included 
in the testing framework developed in Section 2. As a special case, our theory 
applies to the most popular multivariate linear time series models such as 
vector- valued ARMA and linear processes (see Section 4.1). Our focus, how- 
ever, is more on multivariate nonlinear financial time series. As examples, 
we introduce here Bollerslev's [12] constant conditional correlation GARCH 
model (Section 4.2), an extension of this model introduced by Jeantheau [34] 
(Section 4.3), factor GARCH models originally considered in Engle, Ng and 
Rothschild [24] (Section 4.4) and multivariate exponential GARCH models 
(Section 4.5). 

4.1. Multivariate ARMA and linear processes. The theory of linear mul- 
tivariate time series is closely connected to investigating the properties of 
linear processes. Let (ej : j € Z) be a sequence of d-dimensional independent 
and identically distributed random variables with mean and covariance 
matrix ^f. A d-dimensional linear process (yj : j G Z) can then be constructed 
by letting 

oo 

(4.i) yj = J^CtEj-t, j£Z, 

1=0 
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where (Cg :£>0) is a sequence of d x d matrices with absolutely summable 
components. It is clear that (y^ : j S Z) is a strictly stationary but dependent 
process that can be included in the setting of Section 2. 

Theorem 4.1. A multivariate linear process (yj-j € Z) specified by 
(4-1) satisfies Assumption 2.1 if E\\eQ\ A \ < oo and 

E E \ct\«x>. 

m>li>m+l 

The proof of Theorem 4.1 is given in Appendix B. Arguably, the most 
useful class of multivariate linear processes is given by ARMA(j>, q) models. 
These are defined in terms of the linear difference equations 

v q 
(4.2) yj - A-iYi-l = £j + E B t e j-t> 3 e ^ 

where A\,...,A p and B\, . . . , B q are dx d matrices. The following corollary 
identifies those multivariate ARMA models that are special cases of the 
linear process in (4.1). 

Corollary 4.1. A multivariate ARMA(p, q) process (yj :j £ Z) speci- 
fied by (4-2) satisfies Assumption 2.1 if E\\£q\ a ] < oo and det^4(z) ^ for all 

z £ C such that \z\ < 1, where A(z) = 1 — A\z A p z p is a matrix-valued 

polynomial. 

Proof. It follows from Theorem 11.3.1 in [16] that the condition 
det A(z) ^ implies the existence of a unique strictly stationary solution 
to the difference equations (4.2) which is given in the form (4.1). Moreover, 
the {Cg :£ > 0) in (4.1) are uniquely determined by the matrix-valued power 
series 

oo 

C(z)=Y,C e z e = A-\z)B(z), \z\<l, 

1=0 

and decay at a geometric rate, whence we have verified the summability 
condition in Theorem 4.1 and complete the proof. □ 

4.2. The constant conditional correlation GARCH model. In this sub- 
section, we discuss one possible way of specifying multivariate models that 
exhibit heteroscedasticity. It is a nontrivial task to extend the univariate 
ARCH and GARCH settings as introduced by Engle [23] and Bollerslev 
[11] to the vector case. Bollerslev [12] suggested the following constant 
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conditional correlation (CCC) model. Let (sj-.j £ Z) be a sequence of d- 
dimensional independent and identically distributed random vectors with 
mean and covariance matrix ^f. Let the process (yj :j £ Z) be defined by 

(4.3) y j = a j oe j , 

v q 

(4.4) CTj o <Tj = u> + ^ ace <Tj-i ° CTj-t + ^2(3 e o yj ^ e o yj_ e , 

e=i e=i 

where u is coordinatewise strictly positive, a\, . . . ,a p and /3 1 ,...,/3_ are 
coordinatewise nonnegative d-dimensional vectors, and p and q are nonneg- 
ative integers. Moreover, o denotes the Hadamard product of two identically 
sized vectors (or, matrices), which is computed by elementwise multiplica- 
tion. A process given by (4.3) and (4.4) formally resembles the structure of 
a univariate GARCH(p, q) time series. In fact, each coordinate is specified 
by a one-dimensional GARCH equation, whose orders are at most (p,q). 
(Lower orders can easily be achieved through zero coefficients.) 

Theorem 4.2. A multivariate CCC-GARCH process (y^jeZ) speci- 
fied by (4-3) and (4-4) satisfies Assumption 2.1 if E[\e$\ A ] < oo and 

r 

lc = max V \\a^i + /3^eL||2 < 1, 

where r = max{p,q}, a e = (a tjl , . . .,a^ d ) T and (3 e = (/3 i)U . . .,f3i >d ) T . Addi- 
tional coefficients appearing in the latter display are set to equal zero. 

The proof of Theorem 4.2 is given in Appendix B. The main condition, 
7c < 1, ensures that the process (yj :jeZ) has finite fourth-order moments 
and is easy to verify based on a decomposition given in Appendix B. Nelson 
[42] showed that 7c < 1 is also necessary in the case of scalar GARCH(1, 1) 
models. Note, however, that even though necessary and sufficient conditions 
for the existence of fourth-order moments are established in the literature 
(see He and Terasvirta [28]), it is in general difficult to decide whether a 
given parameterization satisfies these requirements. 

To further demonstrate the applicability of our method to the transformed 
process (xj :j 6 Z) defined in (3.1), we show now that (3.2) can be verified 
for the multivariate CCC-GARCH model even if 5 < 1/2. For simplicity of 
the presentation, we assume for the moment that p = q = 1 and that each 
coordinate process y,- , possesses 45 moments (5 < 1/2). Then, Theorem 3 of 
[42] implies that 

7C = max E[(ai + PiS^) 26 ] < 1 
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is necessary and sufficient for E'fjyo.il 45 ] < °°- Moreover, using the GARCH(1, 1) 
representation obtained in [8] (see also [42] and display (B.4) in Appendix B) 

and the truncated variables in (B.5), one finds that, for alH = 1, . . . ,d, 



\\yo,i-yif\\ls= (\ui\E[\e^ 5 ]E 



(oo n \ • 

E II 
n=m+l k=l / 



2<5-i \ 1/4 



< I H£[|e ,i| 45 ]£ 



e ii(«i+A^fc,i) 

.n=m+l k=l 



2S 



1/4 



oc 



= m 



\E[\e ,i\ 45 ] E (^[(oi + ft4 4 ) M ]) r 



1/4 



n=m+l 



using the independence of the (£ji-j £ Z) and the relation \y/a — Vb\ 2 < 
| a — 6| for all a,b > to obtain the first equality sign, and 25 < 1 for the 
first inequality sign. The preceding proves via (3.2) that (2.5) holds for the 
transformed variables Xj = \yj\ s if 5 < 1/2. To employ the break detection 
procedure of Section 2 consequently requires not even a finite variance- 
covariance structure if dealing with a multivariate CCC-GARCH process. 
A similar argument applies also in the finite second moment case 5 > 1/2. 

We conclude this subsection with the remark that the traditional GARCH 
specification in (4.4) can easily be replaced with other specifications that 
take into account, for example, asymmetries, as long as necessary and suf- 
ficient conditions for the existence of strictly stationary solutions and mo- 
ments are available. One such class of GARCH models has been studied in 
Aue, Berkes and Horvath [5]. 



4.3. Jeantheau's constant conditional correlation GARCH model. An ex- 
tension of the CCC-GARCH model was introduced by Jeantheau [34] by 
replacing the coefficients oci, . . . ,oc p and /3 l5 . . . ,/3 q in (4.4) with entrywise 
nonnegative d x d matrices A±, . . . , A p and B\, . . . , B q . This leads to the 
modified definition 

v 1 

(4.5) o-j o a-j = u> + E A e(<7j-e a j-t) + E B t(yj~e ° yj-t) 

e=i e=i 

for the conditional covariances. Therein, to is coordinatewise strictly positive. 
Letting A% = diag(o^) and = diag(/3^), obviously returns (4.4). In this 
case, the criteria for the existence of a unique, nonanticipative (i.e., future 
independent) strictly stationary solution can be obtained mimicking the 
univariate results provided by Brandt [15] and Bougerol and Picard [13, 14]. 
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The more general case (4.5) does not follow straightforwardly and requires 
a number of modifications. 

We introduce more notation. For a G (0,1), we extend the L a -norm by 
letting here || • || Q = E[\ ■ \ a ]. For a > 1, we continue to work with || • || Q = 
(E[\ ■ \ a ]) 1 / a . Notice that the matrix norm \ -\e defined below display (2.6) 
satisfies the condition \MN\e < |M|^|iV|£; for two d x d matrices M and 
N . For a d x d matrix M, let finally 

IIMI 



E,a = \\\Jyl \E\\a 

be the composition of the matrix norm | • \e and the L a -norm || • || Q . Now, 
we can give a sufficient condition for the existence of a strictly stationary 
solution in Jeantheau's model. 

Theorem 4.3. If E[\e \ 2a ] < oo and 

r 

U,a = W A e + B e E o\\E, a < 1 
e=i 

for some a > 0, then the (4-5) have the unique, nonanticipative strictly sta- 
tionary and ergodic solution 

oo k 

a j oa J =u+ J2 Yl H( A ie+ B k E j-h h) ^> 3 G Z > 

-k=ll<ii<---<ik<r£=l 

where r = m&x{p, q} and Ej = diag(ej o Sj). Additional matrices appearing 
in the definition of 7j iQ are set to equal the zero matrix. 

With this result at hand, it can further be established under which con- 
ditions the testing procedures of Section 2 are applicable for Jeantheau's 
model. The proofs of Theorems 4.3 and 4.4 are given in Appendix B. 

Theorem 4.4. A multivariate CCC-GARCH process (yj-.jeZ) spec- 
ified by (4-3) and (4-5) satisfies Assumption 2.1 if E^eq^] < oo and jj = 

U,2 < 1. 

Observe, however, that 7j < 1 here and 7c* < 1 in Bollerslev's [12] CCC- 
GARCH model are not even equivalent if Ag = di&g(ai) and E>i = diag(/3^) 
are diagonal matrices. Rather, the condition 7c < 1 is less restrictive. This 
can be seen from the inequality 

r r 

>yj = V \\Ai + B £ E \\e,2 > max V \\ae ti + /3^£n ih =7C- 

Standard arguments imply that this inequality is strict if, for example, £0 is 
normally distributed (or, if it has a positive density on M. d ). 
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4.4. Dynamic factor models. It is often believed in economic theory that 
the dynamics of a multivariate sequence (y^ :j£Z) can be adequately de- 
scribed by so-called (unobserved) common factors, say, (zj : j £ Z) with val- 
ues in R rf ' , d' < d. This gives rise to the model 

(4.6) y j = Lz j +Z j , j€Z, 

where L is & d x d' matrix of factor loadings and (£j :j 6 Z) a sequence 
of tf -dimensional errors. The following theorem specifies conditions under 
which the sequence (y^ :j€Z) satisfies Assumption 2.1. 

Theorem 4.5. Let (z,,- : j 6 Z) and (£j :j£Z) 6e suc/i i/iai t/iey satisfy 
Assumption 2.1. Then, the process (y,-:j'GZ) specified by (4-6) satisfies 
Assumption 2.1 as well. 

The main argument for the use of factor models is their potential to reduce 
dimensionality and consequently the number of parameters to be estimated 
from data. Often, the components of the factors are assumed to be uncor- 
rected, indicating genuinely different driving forces that account for the 
overall volatility behavior. The first factor GARCH model was introduced 
by Engle, Ng and Rothschild [24] , more recent factor models were put forth 
in [38, 49]. For an in-depth overview we refer to [47]. 

4.5. Multivariate exponential GARCH. A further subclass of the multi- 
variate GARCH family is provided by Kawakatsu's [35] matrix extension of 
Nelson's [42] exponential GARCH model. The model in [35] is given by first 
assuming that the random vectors (y^ :jGZ) follow the equations 

(4.7) y, "] J'eZ, 

with (ej :j S Z) being a sequence of independent, identically distributed d- 
dimensional random vectors. The d x d matrices Hj are now specified by 
letting 

p 

vech[log Hj -C]=^2Ai vech[log Hj_e - C] 

(4.8) 

q q 
+ B e £ j-e + F d\ £ j-i\c ~ E[\ej-i\ c ]), 
i=i i=i 

where C denotes a symmetric dxd matrix and | • | c component-wise absolute 
value. The parameter matrices A#, Bp and F% are of dimensions d x d, d x d 
and dxd, respectively. Since log Hj is clearly a symmetric matrix, we have 
that 

OO -i 

(4.9) Hj = exp(\ogHj) = £ ^ (log Hj) k , j € Z, 
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is a positive definite matrix. The apparent advantage of this approach is that 
no (unnatural) restrictions on the parameters have to be imposed through 
the model specification so as to ensure positive definiteness. The dynamic 
correlation structure is enabled through the recursions in (4.8) which de- 
termine logHj, while Hj in turn is obtained from the matrix exponential 
representation in (4.9). We shall work here with a slight modification given 
by 

(4.10) vechflogiTj -C}= Avechflogi^i - C] + F(e,-_i, 

where A is a d x d matrix and F is a measurable function taking values in 
R a . The following theorem provides a sufficient condition for the existence 
of a strictly stationary solution to the multivariate exponential GARCH 
equations. Denote by log + x = maxjlogx, 0}. 

Theorem 4.6. If E[log + \F(e q , . . . ,ei)|] < oo and 

(4.11) limsupilog|A fc | s < 0, 

k— >oo ™ 

then the (4-10) have the unique, nonanticipative strictly stationary and er- 
godic solution 

oo 

vechflogFj - C] = £ A k F(ej- k -i, . . .,ej- k - q ), j £ Z. 

fc=0 

Condition (4-H) is in particular satisfied if \A\e < 1. 

To prove that Assumption 2.1 is applicable here, we need to assume ad- 
ditionally that the random variable \F(e—i, . . . ,S—q)\ has a finite moment 
generating function in a sufficiently large neighborhood of zero. Details are 
provided next. 

Theorem 4.7. A multivariate exponential GARCH process (y,-:_7'€Z) 
specified by (4-7) and (4-10) satisfies Assumption 2.1 if \A\e < 1 and 
-E[exp(t|F(e_i, . . . , e_ g )|)] < oo for some t > y/8q. 

The proofs of Theorems 4.6 and 4.7 are postponed to Section B.4. 

4.6. Summary. In this section, we have discussed the most popular mul- 
tivariate GARCH models. Which of these models a practitioner will even- 
tually pick will depend heavily on the application at hand. All specifica- 
tions have their particular strengths and weaknesses. Research — applied and 
theoretic — is still to be conducted to reveal further insight. We refrain hence 
from giving recommendations here but refer to the survey papers [7, 47] for 
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Table 1 

Critical values for Q(X>) for the three significance levels 0.10, 0.05 and 0.01. For each 
critical value, we report P(f2(0) < zi- a ) in brackets. The last column, = oo, 
corresponds to the standard normal case as stated in Remark 2.1 






10 


15 


20 


50 


100 


200 


500 


oo 


aa9o(0) 


1.33 


1.33 


1.32 


1.31 


1.31 


1.30 


1.29 


1.28 




(0.89) 


(0.89) 


(0.89) 


(0.89) 


(0.90) 


(0.90) 


(0.90) 


(0.90) 


«6W*0 


1.84 


1.81 


1.79 


1.74 


1.71 


1.69 


1.68 


1.64 




(0.93) 


(0.94) 


(0.94) 


(0.94) 


(0.94) 


(0.94) 


(0.95) 


(0.95) 


<?0 n 99(0) 


2.90 


2.80 


2.74 


2.59 


2.51 


2.46 


2.41 


2.33 




(0.98) 


(0.98) 


(0.98) 


(0.98) 


(0.99) 


(0.99) 


(0.99) 


(0.99) 



additional reading. Independent of practical modeling considerations, the 
tests introduced in Section 2, however, are flexible enough to cover any of 
the particular multivariate specifications given in Sections 4.1-4.5. To inves- 
tigate their finite sample behavior will be subject of Section 5. 

5. Empirical results and applications. The empirical part of the paper 
includes in its first part a discussion on the computation of the critical values 
and the validity of the normal approximation provided in Remark 2.1 in the 
case of finite sample sizes (Section 5.1). In a second part, we assess the finite 
sample properties of the test procedures for several multivariate time series 
models (Section 5.2). Finally, we address the applicability by investigating 
volatilities and cross- volatilities of stock data for 12 companies (Section 5.3). 

5.1. Computation of critical values. The asymptotic critical values of 
the test statistics A n and Q, n used in Theorem 2.1 can be computed either 
using the precise distribution of the limits A(t)) and 0(D), respectively, or — 
provided X) is sufficiently large — the normal approximation given in Remark 
2.1. Critical values for d = 1,...,5 have already been tabulated in Kiefer 
[36]. Here, we provide critical values for larger values of t) that seem to be 
more realistic in today's financial applications. The calculations are based 
on the limit distributions of A(i)) and Q(d) which have been derived in [36]. 

Let <7p_a,(t>) be the 1 — a quantile of 0(t>)> the standardized version of 
0(c)), and let furthermore qi , _ a (oo) = z\- a denote the 1 — a quantile of the 
standard normal distribution. In Table 1, we have recorded the critical values 
for three significance levels, namely a = 0.10, 0.05 and 0.01, and various 
values of d ranging from 10 to 500. Under each critical value, we report in 
brackets also the values P(£l(d) < zi- a ). Let qi- a (P) be the 1 — a quantile 
of A(D), the standardized version of A(fl). The corresponding critical values 
for the same choices of a and d have been collected in Table 2. 
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It can be seen from Table 1 that the normal approximation (0 = oo) works 
very reasonable in the case of fi(fl). Even for those quantiles q^_ a that are 
the farthest from the corresponding standard normal quantile z±- a , it still 
holds that P(Q(d) < Zi- a ) wl-a. This, however, does not hold true for 
the second statistic A(t)). Here, the normal approximation (U = oo) works 
considerably worse and it cannot even be recommended for = 500. The 
exact values should consequently be used instead. 

5.2. A simulation study. In this subsection, we report the results of a 
small simulation study that includes multivariate AR(1), dynamic factor 
model and exponential GARCH specifications for the process (y~j :j£Z). In 
all cases, we use the statistics Q n . Due to the weak dependence displayed 
by the data generating processes, we choose to work with the Bartlett es- 
timator as a proxy for the asymptotic covariance matrix X with window 
length q = q(n) = log 10 n (see [4]). It is shown in Chapter 11 of [16] that the 
Bartlett estimator satisfies condition (2.6). 

If ( v j : 3 ' G ^) follows a parametric model, then X is a function of the model 
parameters. These, in turn, yield a natural plug-in estimator X* based on 
the relevant parameter estimators. This method, however, fails to work if the 
model is not correctly specified. Then, the plug-in estimator will converge 
to a limit X* ^ X, resulting in a potential power loss of the test statistics. 
Moreover, estimation procedures and goodness of fit tests for the nonlinear 
models discussed in Section 4 still often lack accuracy and a nonparametric 
approach seems — at this point — more advisable. 

Autoregressive PROCESSES. We used as data generating process the 
four-dimensional AR(1) time series given by the equations 

yj =A yj-l + £ ji 3 = l,---,n, 

Table 2 

Critical values for A(0) for the three significance levels 0.10, 0.05 and 0.01. For each 
critical value, we report P(A(3) < zi- a ) m brackets. The last column, = oo, 
corresponds to the standard normal case as stated in Remark 2.1 



D 


10 


15 


20 


50 


100 


200 


500 


oo 


9a 9 o(0) 


2.64 


2.53 


2.46 


2.27 


2.16 


2.06 


1.96 


1.28 




(0.56) 


(0.58) 


(0.59) 


(0.63) 


(0.66) 


(0.69) 


(0.72) 


(0.90) 


?o. 9S (0) 


3.17 


3.02 


2.92 


2.69 


2.55 


2.44 


2.33 


1.64 




(0.69) 


(0.71) 


(0.72) 


(0.76) 


(0.78) 


(0.81) 


(0.83) 


(0.95) 


?0. 99 (f) 


4.28 


4.04 


3.89 


3.53 


3.33 


3.18 


3.04 


2.33 




(0.85) 


(0.87) 


(0.88) 


(0.91) 


(0.93) 


(0.94) 


(0.95) 


(0.99) 
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where yo has been obtained after a burn-in phase of 500 iterations. The 
sample sizes under consideration are n = 200, 500, 800 and 1000. The in- 
novations £i,...,e n have been specified as independent 4-variate normal 
variates with independent components. To assess the empirical level of the 
test, we have chosen A = (0.1)/, where I denotes the 4x4 identity matrix. 
Each component of the process is therefore a univariate AR(1) time series 
which is independent of the others. 

To assess the power, we introduce several alternatives. These are specified 
by the form of the parameter matrix before and after the change-point k* = 
n/2. In particular, we work here with A = (0.1)/ prior to k* and with A* = 
(0.1)7 + 51 after k*, where 1 denotes the 4x4 matrix for which all entries 
are equal to 1, thus introducing correlation between the coordinates after 
k* . The values of 5 vary between 0.1 and 0.6. Using Theorem 2.1, the critical 
values have been computed using the exact distribution of the limit random 
variable f2(Z)) given in (2.9) and discussed in the previous subsection. All 
results are based on N = 1000 repetitions. 

The results are summarized in Table 3 for the three levels a = 0.10, 0.05 
and 0.01. The row with 5 = corresponds to the empirical levels of the 
test. It can be seen that the procedure is conservative in finite samples 
as the empirical levels stay below the asymptotic critical values. Empirical 
and asymptotic level are close for the sample size n = 1000. On the other 
hand, small changes are found with greater difficulty if the sample size is 
small (n = 200). The power for bigger changes and/or larger sample sizes is, 
however, very reasonable. 

Factor models. As second data generating process we used a factor 
model in which the four-dimensional yj are given by the equations 

y j = Lz j + £ j , j = l,...,n, 

Table 3 

The rejection levels of the statistic Q n for the various AR(1) processes of Section 5.2. 
The row specified by 5 = contains the empirical levels 



a 0.10 0.05 0.01 



S/n 


200 


500 


800 


1000 


200 


500 


800 


1000 


200 


500 


800 


1000 


0.0 


0.07 


0.09 


0.10 


0.10 


0.03 


0.04 


0.04 


0.04 


0.00 


0.01 


0.01 


0.01 


0.2 


0.22 


0.42 


0.67 


0.78 


0.11 


0.29 


0.53 


0.66 


0.03 


0.11 


0.25 


0.36 


0.3 


0.40 


0.76 


0.94 


0.98 


0.25 


0.61 


0.89 


0.96 


0.07 


0.31 


0.72 


0.86 


0.4 


0.71 


0.96 


1.00 


TOO 


0.57 


0.91 


1.00 


1.00 


0.26 


0.74 


0.98 


1.00 


0.5 


0.94 


1.00 


1.00 


TOO 


0.90 


1.00 


1.00 


1.00 


0.73 


0.97 


1.00 


1.00 


0.6 


1.00 


1.00 


1.00 


TOO 


0.99 


1.00 


1.00 


1.00 


0.98 


TOO 


1.00 


1.00 
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where L is a 4 x 2 matrix, £i> • ■ • >£n are independent zero mean 4-variate 
normal variates with unit covariance matrix. The 2-variate factors zi, . . . ,z n 
are given by the CCC-GARCH equations (4.3) and (4.4) with specifications 

fl.0\ /0.3\ /0.3 



i.oy ' 1 V - 3 /' V - 3 

and bivariate normal Sj with identity covariance matrix. To assess the em- 
pirical level, we chose 



L 



(I 0\ 

1 

1 

Vo \J 



under Hq, while for the power considerations we picked L* = 5L under Ha- 
The choices for the number of repetitions N, the sample sizes n and the 
significance level a are the same as for the AR(1) processes above, while for 
the change sizes we used values of 5 ranging from 1.1 to 2.0. The results are 
summarized in Table 4. The empirical levels appear to be more volatile than 
in the AR(1) case but quite satisfactory for the sample size n = 1000. The 
findings for the power of the test are similar to what has been observed for 
the AR(1) time series. 

Exponential GARCH models. The final data generating process is 
a four-dimensional exponential GARCH model given by the (4.7) and (4.10) 
with specifications p = q = 1, A\ = 0.1 diag(I( 10 )), B T = (0.1 diag(I( 4 )), 0) 
and F = 0, where \rj\ denotes the d-dimensional vector whose entries are 
all equal to 1 and the 4x4 zero matrix. Under Hq, we worked with 
C = 0.2diag(I( 4 )), while under Ha we used C prior to k* and with C* = 5C 
thereafter, where 5 ranges from 1 to 3.5. The errors {sj'-j = l,...,n) are 

Table 4 

The rejection levels of the statistic Q n for the factor model of Section 5.2. The row 
specified by 5 — 1 contains the empirical levels 

a 0.10 0.05 0.01 



S/n 


200 


500 


800 


1000 


200 


500 


800 


1000 


200 


500 


800 


1000 


1.0 


0.15 


0.12 


0.13 


0.15 


0.06 


0.05 


0.07 


0.08 


0.01 


0.01 


0.01 


0.01 


1.1 


0.16 


0.19 


0.26 


0.31 


0.09 


0.09 


0.16 


0.21 


0.02 


0.02 


0.05 


0.06 


1.2 


0.25 


0.39 


0.58 


0.66 


0.13 


0.27 


0.45 


0.55 


0.03 


0.12 


0.23 


0.30 


1.4 


0.52 


0.85 


0.96 


0.99 


0.37 


0.75 


0.93 


0.98 


0.13 


0.51 


0.83 


0.92 


1.6 


0.76 


0.98 


1.00 


1.00 


0.63 


0.96 


1.00 


1.00 


0.35 


0.86 


0.99 


1.00 


1.8 


0.91 


1.00 


1.00 


1.00 


0.82 


1.00 


1.00 


1.00 


0.56 


0.99 


1.00 


1.00 


2.0 


0.96 


1.00 


1.00 


1.00 


0.93 


1.00 


1.00 


1.00 


0.74 


0.99 


1.00 


1.00 
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Table 5 

The rejection levels of the statistic fl n for the exponential GARCH model of Section 5.2. 
The row specified by 5 = 1 contains the empirical levels 



a 0.10 0.05 0.01 



S/n 


200 


500 


800 


1000 


200 


500 


800 


1000 


200 


500 


800 


1000 


1.0 


0.07 


0.07 


0.07 


0.08 


0.04 


0.03 


0.04 


0.04 


0.00 


0.01 


0.00 


0.01 


1.5 


0.10 


0.19 


0.26 


0.30 


0.04 


0.09 


0.16 


0.18 


0.00 


0.01 


0.04 


0.06 


2.0 


0.18 


0.47 


0.75 


0.83 


0.08 


0.31 


0.61 


0.74 


0.01 


0.11 


0.34 


0.49 


2.5 


0.40 


0.83 


0.98 


0.99 


0.25 


0.74 


0.96 


0.98 


0.07 


0.47 


0.85 


0.95 


3.0 


0.61 


0.98 


1.00 


1.00 


0.43 


0.96 


1.00 


1.00 


0.17 


0.86 


1.00 


1.00 


3.5 


0.83 


1.00 


1.00 


1.00 


0.70 


1.00 


1.00 


1.00 


0.38 


0.98 


1.00 


1.00 



standard normal and the choices for n, k* and a are the same as before. 
The results are summarized in Table 5. They are close to the findings in the 
AR(1) and factor model cases. 

Estimating the break location. To conclude the simulation study, 
we briefly investigated also how well the relative change location is ap- 
proximated by the estimator n of (2.12) for the AR(1) processes, the factor 
and exponential GARCH models. In Table 6, we have reported mean, me- 
dian and standard deviation for the various choices of 5 and n specified in 
the foregoing paragraphs and for the break point k* = n/2, thus leading to 
= 0.5. The significance level was set to a = 0.1. The table shows that for all 
three processes there is a positive bias which diminishes as the sample size 
or the change magnitude 5 increase. Generally, using the median leads to a 
superior approximation of the true than utilizing the mean, thereby indi- 
cating the existence of unusually large deviations (outliers) in several of the 
simulation runs. The estimator n works satisfactory in the case n = 1000. 
This shows that the asymptotic result in display (2.13) applies reasonably 
well here, although future research should pursue a more detailed in-depth 
analysis of the performance of n . 

5.3. An application to financial data. We study in our application the log 
returns of the adjusted closing stock prices between July 19, 1993, and March 
19, 2009, of the 12 companies listed in Table 7, yielding the 12-dimensional 
observations yi, • • • ,y394i- There are four companies in the airline sector, 
four in the automotive sector and four in the energy sector. More precisely, 
if pj t £ denotes the price of stock £ at time j, then we work with the centered 
log-returns 

/ \ 1 3941 , \ 

«'' = l0g U7)-394lg l0g U7j' 
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Table 6 

The estimation of the change-point fraction 6 for the AR(1) processes (upper panel), the 
factor GARCH model (middle panel) and the exponential GARCH model (lower panel) 



n 
S 




200 






500 






800 






1000 




mean 


med 


sd 


mean 


med 


sd 


mean 


med 


sd 


mean 


med 


sd 


0.2 


0.55 


0.55 


0.11 


0.56 


0.55 


0.09 


0.54 


0.53 


0.08 


0.54 


0.52 


0.07 


0.3 


0.56 


0.55 


0.10 


0.55 


0.53 


0.08 


0.54 


0.52 


0.06 


0.53 


0.52 


0.05 


0.4 


0.59 


0.57 


0.09 


0.55 


0.53 


0.06 


0.53 


0.52 


0.05 


0.53 


0.51 


0.04 


0.5 


0.58 


0.57 


0.08 


0.54 


0.52 


0.05 


0.53 


0.52 


0.05 


0.53 


0.51 


0.04 


1.2 


0.56 


0.56 


0.11 


0.54 


0.53 


0.09 


0.54 


0.52 


0.09 


0.53 


0.52 


0.08 


1.4 


0.56 


0.55 


0.09 


0.53 


0.52 


0.07 


0.52 


0.51 


0.05 


0.52 


0.51 


0.04 


1.6 


0.55 


0.53 


0.07 


0.52 


0.51 


0.05 


0.51 


0.51 


0.03 


0.51 


0.51 


0.03 


1.8 


0.53 


0.52 


0.06 


0.52 


0.51 


0.03 


0.51 


0.51 


0.02 


0.51 


0.50 


0.02 


2.0 


0.53 


0.52 


0.05 


0.51 


0.50 


0.03 


0.51 


0.50 


0.02 


0.51 


0.50 


0.02 


1.5 


0.52 


0.54 


0.12 


0.54 


0.53 


0.09 


0.52 


0.52 


0.09 


0.53 


0.52 


0.09 


2.0 


0.54 


0.53 


0.10 


0.53 


0.52 


0.07 


0.52 


0.51 


0.07 


0.52 


0.51 


0.06 


2.5 


0.54 


0.53 


0.09 


0.52 


0.51 


0.04 


0.51 


0.50 


0.04 


0.51 


0.50 


0.03 


3.0 


0.53 


0.52 


0.07 


0.51 


0.50 


0.02 


0.51 


0.50 


0.02 


0.51 


0.50 


0.02 


3.5 


0.53 


0.51 


0.06 


0.51 


0.50 


0.02 


0.51 


0.50 


0.02 


0.50 


0.50 


0.01 



j = l,...,3941;£ = l,...,12. 

To illustrate the time- varying nature of the volatilities, we have computed 
the rolling volatility and cross-volatility estimators 

1 j 

x fj (k,£) = — y^Vi,i, j = 101,..., 3941; M = l,... ,12. 

iUU i=j-100+l 

For each of the three sectors, the corresponding estimators are given in Fig- 
ure 1. They appear to be time-dependent. To further assess this conjecture, 
we computed the test statistic value ^3941 = 60.07. Since d = 12, we have 
d = 78. This leads to the critical value 12.00 at the 95% level. Therefore, 
the hypothesis of no change in the volatility and cross- volatility structure is 
rejected. 

To illustrate the potential to detect multiple departures of the stability 
hypothesis, we have applied a binary segmentation procedure of the original 
test. Each time we reject Hq, we re-apply it to the subsamples obtained 
from splitting the data into two pieces at the observation determined by 
(2.12), that is, k* = \ 6 n n\ . The corresponding findings are reported in Table 
8. A number of the detected changes can be readily associated with major 
historical events. For example, the change found in the first round of the 
segmentation procedure along with the other estimated changes in 2001 are 
linked to the bursting of the dot-com bubble and the September 11 attacks, 
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1994 1996 1998 2000 2002 2004 2006 2008 




1994 1996 1998 2000 2002 2004 2006 2008 




1994 1996 1998 2000 2002 2004 2006 2008 1994 1996 1998 2000 2002 2004 2006 2008 



q - 




1994 1996 1998 2000 2002 2004 2006 2008 1994 1996 1998 2000 2002 2004 2006 2008 



Fig. 1. The volatilities (left) and cross-volatilities (right) of the airline sector (upper 
panel), the automotive sector (middle panel), and the energy sector (lower panel). The 
company abbreviations used in the plots are detailed in Table 7. 
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Table 7 

The 12 companies whose stock values are studied in Section 5.3 





Symbol 


Name 


Sector 


1 


ALK 


Alaska Air Group, Inc. 


Airline 


2 


AMR 


AMR Corp. (American Airlines) 


Airline 


3 


CAL 


Continental Airlines, Inc. 


Airline 


4 


LUV 


Southwest Airlines Co. 


Airline 


5 


F 


Ford Motor Co. 


Automotive 


6 


GM 


General Motors Corp. 


Automotive 


7 


HMC 


Honda Motor Co. 


Automotive 


8 


TM 


Toyota Motor Corp. 


Automotive 


9 


APA 


Apache Corp. 


Energy 


10 


APC 


Anadarko Petroleum Corp. 


Energy 


11 


OXY 


Occidental Petroleum Corp. 


Energy 


12 


RIG 


Transocean Inc. 


Energy 



while the break dates in 1997 and 1998 are connected to the Asian financial 
crisis, and the collapse of the hedge fund Long-Term Capital Management 
and the Russian financial crisis, respectively. The detected breaks in 2007 
and 2008 can be related to the collapse of the housing market in the United 
States and several European countries. Finally, the latest significant break, 
which was found at lag 3809, signifies the beginning of the recent financial 

Table 8 

The segmentation procedure performed on the entire data set. The estimated 
change-points are listed with their corresponding test statistic value, the round in which 
the change was found and whether or not it was significant at the 95% error level 



k* 


k* (Date) 




Round 


Significant 


305 


1994-09-30 


11.60 


4 


No 


647 


1996-02-07 


19.80 


3 


Yes 


1021 


1997-07-31 


14.78 


4 


Yes 


1259 


1998-07-13 


34.05 


2 


Yes 


1400 


1999-02-02 


13.60 


4 


Yes 


1770 


2000-07-25 


18.65 


3 


Yes 


1886 


2001-01-11 


12.07 


4 


Yes 


2101 


2001-11-23 


60.07 


1 


Yes 


2358 


2002-12-03 


14.66 


4 


Yes 


2728 


2004-05-24 


27.90 


3 


Yes 


3175 


2006-03-03 


14.70 


4 


Yes 


3589 


2007-10-24 


28.95 


2 


Yes 


3710 


2008-04-18 


10.11 


4 


No 


3809 


2008-09-09 


12.13 


3 


Yes 


3871 


2008-12-05 


7.47 


4 


Xo 
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crisis. The detected change-point time is September 9, 2008, thus predating 
the collapse of the investment bank Lehman Brothers by only three trading 
days. 

APPENDIX A: A MULTIVARIATE FUNCTIONAL 
CENTRAL LIMIT THEOREM 

In this section, we derive a general multivariate functional central limit 
theorem (FCLT) and state moreover a derivative that is tailor-made for 
the examples discussed in Section 4. The results extend the corresponding 
univariate theory provided, for example, in Billingsley's monograph [10]. 

Recall from Assumption 2.1 that (sj :j £ Z) is a sequence of independent, 
identically distributed random variables with values in M d and that we fur- 
ther suppose that the (/-dimensional sequence (y,- :j G Z) is given by the 
relations (2.3). To prove a FCLT for this sequence, we shall utilize the se- 
quence of m-dependent random vectors (y^ :jGZ) given in (2.4), which is 
sufficiently close to (yj £ Z) in the sense of condition (A.l) below. Recall 
finally that a stochastic process (Wp(t) :t £ [0, 1]) is called a (/-dimensional 
Brownian motion with covariance matrix T £ M. dxd if it is Gaussian with 
mean and covariance function Cov(Wr(s), Wp(i)) = min{s,t}r. 

Theorem A.l. Assume that the d- dimensional random process (yj :j € 
Z) specified by (2.3) satisfies E[yj] =0 and E[\yj\ 2 ] < oo. Suppose further 

that, for any m>l, the vectors y^ can be defined such that 

(A.i) £llyo-yi m) ll 2 <~- 

m>l 

Then the series T = J2jez Cov(yo, yj) converges (coordinatewise) absolutely 
and 

i [nt] 

1 ^-i D d [0,l] . s , . 

vn j=i 

where the convergence takes place in the d- dimensional Skorohod space Dd[0, 1]. 
The proof of Theorem A.l needs the following small auxiliary result. 

Lemma A.l. Let (X n :n> 1) and (Y n :n> 1) be sequences of random 

vectors in W 1 such that X n — > X and Y n Y as n — > oo, where the 
limit random vectors X and Y are independent. Suppose further that there 
are sequences (X* :n > 1) and (Y* :n > 1) such that (a) X* and Y* are 

independent for all n > 1, and (b) |X* — X n | and |Y* — Y n | as 

n— >oo. Then it holds that (X. n ,Y n ) (X.,Y) asn^oo. 
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Proof. By assumption, it holds that X* — ► X and Y* — ► Y. Since 
the sequences (X* :n > 1) and (Y*:n > 1) are independent, we obtain 

(X*,Y*) (X,Y). Viewing (X*,Y*) and (X,Y) as random vectors in 

M. 2d implies | (X* , Y*) — (X n , Y n )| — > and consequently the assertion. □ 



Proof of Theorem A.l. 

Step 1. For n > 1 and t G [0,1], let V n (t) = n~ 1 / 2 E^iYj- It will be 
shown that the finite-dimensional distributions of the partial sum process 
(Vn(t) :t G [0, 1]) converge to the finite-dimensional distributions correspond- 
ing to the d-dimensional Brownian motion (Ws(t) :t G [0, 1]). Let K > 1 and 
< to < t\ < ■ ■ ■ < tx < 1- Because the limit process has independent incre- 
ments, proving the latter convergence is equivalent to showing that 

(V n (t e -i,t e ):l<£<K) 

(A.2) 

-A(Wrfe)-W r (t w ):l<K^) (rwoo), 

where V n (s,t) = n _1//2 E^Lj+i Yj for < s < i < 1. To do so, we focus in 
the first part of the proof on one element of the random vectors in (A.2). 
More precisely, we shall establish with the Cramer- Wold device that, for 
any < s < t < 1, 

(A.3) V„(a 1 i)-^W r (t)-Wr(s), (rwoo). 

For a = (a±, . . . , ad) T G K d , let the univariate sequences (zj :jeZ) and -j G 
Z) be defined by letting 



a^ yj and z y -" 1 ' = a T y^ 



j 

respectively. Let y = (y ,i, • ■ -,yo,d) T and VcT' = G/cu>- ■ -^04 ) T and ob " 
!/o,t - I - |y° _ ^o m) I for all 1 < i < d. Tl 

En ( m )n i in ( m )n 

lko-«o II2 < 2^ maxKlllyo-Yo II2 

m>l — — 
\i=l / V m>l 



serve that |yo,i — — |yo ~~ ycT^I f° r all 1 < £ < d. Therefore, 

y ( m )\\ <■ \ " in,, „M| 

'0 lb - / t - 
m>l m>l 

( m ) 1 1 

12 <°° 



by assumption (A.l). Theorem 21.1 in [10] consequently yields that 

[nt] 

-= E ^-^A(0,a 2 2 (t- S )) (n-oo), 
v n j= [„ s]+ i 
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where a 2 z = J2jez E[ z o z j] = &T Cov(yo, Yj)a < oo, which implies in par- 
ticular that r converges coordinatewise absolutely. On the other hand, it 
holds that 

a T (W r (i) - W r (s)) = N(0, a 2 z (t - s)), 

so that (A. 3) is proved. However, since the partial sums V n (i^_i,fy), 1 < 
£ < K, are dependent, (A. 3) does not yet imply (A. 2). In the second part 
of the proof, we shall therefore construct independent random variables to 
replace them. For < s < t < 1, let V*(s,t) = ra^ 2 £^ ] M+1 yj M+1_i) , 

where y^ = 0. Then V*(s,t) and V n (0, s) are independent. Moreover, for 
any 5 > 0, it holds that 

P(|V„(M) - K(M)I > S) < ±E{\V n (s,t) - V*(s,t)l 2 ] 




\m=0 / 

as n — > oo. An application of Lemma A.l shows that the vectors (V n (0, s),V n (s, 
t)) converge in distribution to (Wp(s), Wr(t) - Wp(s)) as n — > oo. This 
proves (A. 2) for K = 2. The same arguments apply also for arbitrary K. 

Step 2. It remains to verify that the partial sum process (V n (t) : t € [0, 1]) 
is tight. Since it suffices to show the tightness coordinatewise, a reference to 
Theorem 21.1 in [10] completes the proof of Theorem A.l. □ 

The main results in Section 2 cover the case of detecting breaks in the 
covariance structure of multivariate financial time series. Therefore, at least 
in the context of the present paper, the following derivative of Theorem A.l 
is of greater importance. 

Theorem A. 2. Suppose that the assumptions of Theorem A.l hold true 
and that E[\yj\ 4 ] < oo. If (2.5) is satisfied, then the series E = 
J2jez Cov(vech[yoyol, v ech[yjyJ]) converges (coordinatewise) absolutely and 

1 M 

vn j=i 

where d = d{d+ l)/2 and the convergence takes place in the D-dimensional 
Skorohod space D$ [0, 1] . 
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Proof. First, we define the symmetric matrix Yq = yoYo — ^[yoYo ] • ^° 
prove the assertion of Theorem A. 2, it suffices to show that the assumptions 
of Theorem A.l are satisfied for 

z = vech[Y ]- 

Note that the moment assumptions E[zq] = and -E[|zo| 2 ] < oo are clearly 
fulfilled, so that it remains to verify condition (A.l). Recall that yo = 

(yo,l,---,yo,d) T and y m) = (yfo ,y$$) T . The (kj)th element of Y is 
thus given by 

*o(M) =yo,kyo,e - E[y 0jk y 0j i], l<k,£<d. 

In a similar fashion, we define the matrix YJ, corresponding to the approx- 
imating variable Zq™^ = vech[Yg m ^] by letting 

F M (M) = y^y^f - E[y oM ], l<k,£<d. 

Observe next that, applying Cauchy-Schwarz and Minkowski, with some 
constant c > 0, 

j^n fm) (m)|2i 

E[\yo,kyo,£-y 0jk y /\ \ 

, r-,r/i ("'Oil I , I (m)iN2 n 

< ^L(|l/o,fc - I/o,fc M3/o^| + l2/o,* lll/o,* -J/o/l) J 

, rn 1. 2 ii (m)||2-|/|| (m)\l i II ( m )ll \2 

- 1™ £<^" yo ' £ " 4 ' " y °' fc H4}(llyo,fc - %,fc IU + \\vo4 - Vo/ u) 
- 4 1 3 na 5 illfo^lli (llyo.fclk + ho,k - i/o!?IL) 2 }llyo - yo m) ll 2 

l<kl<d 



^ A II II I V II \ m > II II l m i II 2 

< 4 HyolU + 2^ llyo ~ yo u) llyo — yo u 



2 

(m)u \ ii (m)||2 



2 n ( m )l|2 

c ||yo - yo IU' 



II ( m )ll ^ II ("Oil ^ V^OO || ("Oil ^ rpi 

since ||y ,fc - 2/o,fc Ik < llyo ~ Yo Ik < L m =il|yo - y Ik < oo. Thus, we 
arrive at 

oo / d d \ V 2 

Z ll z o - z o m) ll 2 < Z ZZ^[|yo,fcyo^-yo?2/o?l 2 ] 

m>l m=l\fc=l£=l / 



< cd Zlly°~yo m) |l4< 



:3C 



m=l 



and the proof is therefore complete. □ 
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APPENDIX B: MATHEMATICAL PROOFS 
B.l. Proofs of Theorems 2.1, 2.2 and Remark 2.1. 



Proof of Theorem 2.1. Theorem A. 2 yields that 

-= £ (vech[y jy J] - E[vech[y jy J]]) D ^ 1] W s (i) (n - oo). 



Hence, 

Vn Vj=i i=i / 

where (Bs(i):i G [0,1]) denotes a Gaussian process with mean function 
£?[Bs(t)] = and covariance function E[Bz(t)~B'^(s)] = E(min{t, s} — ts). 
Now Theorem 2.1 follows from the continuous mapping theorem and as- 
sumption (2.6). □ 

Proof of Remark 2.1. We consider A(Z>) first and show in a first step 
that 

(B.l) A;(i) = i= 5>f (i) - t[l - t]) °H ] A*(t) (H - oo), 

where C ^iJ indicates weak convergence in the space of continuous functions 
defined on the interval [0, 1] and (A*(t) : t £ [0, 1]) denotes a continuous Gaus- 
sian process. Note that in the discourse of the proof, we shall not need the 
covariance structure of this process and we consequently do not compute it. 
We are going to apply the Cramer-Wold device. It follows from the inde- 
pendence of the Brownian bridge processes (Bg(t) : t £ [0, 1]), 1 < £ < d, that, 
for all U £ [0, 1] and Aj £ E, l<i<N, 

N 

5>a;(^)^£ (a-oo), 

i=l 

where £ denotes a centered normal random variable. Next, we verify the 
validity of the moment condition (12.51) in Billingsley [10]. Let fJ,(t,s) = 
t(l — t) — s(l — s) and denote in the following by c a universal constant 
which may vary from line to line. Observe that Rosenthal's inequality (see, 
e.g., page 59 of Petrov [44]) implies that 

E[(A;(t) - a;(,)) 4 ] 

<"\j2E[(Bnt)-Bl(s)-^t,s)f] 
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k£=1 



<c(t-s) 2 , 



since, by the Cauchy-Schwarz inequality and the fact that E[(Be(t) — Be(s)) 8 ] < 
c(t — s) 4 , we have after routine computations that 

E[(B}{t)-B}{s)-^t,s)) A ]<c{t-s)\ 

E[{B}{t)-B}{s)-n{t,s)f\<c\t-s\. 

Theorem 12.3 in [10] yields now the tightness and the continuity of the limit 
process (A*(t) : t 6 [0,1]). 

Let £j, denote the location of the (first) maximum of B 2 (t) + • • • + B 2 (t). 
Since the function t \— > t(l — t) reaches its largest value for t = 1/2, the weak 
convergence in (B.l) implies that 

11 -0 P (d~ 1/4 ) (B->oo). 



(B.2) \t, 2 

We show now that 



(B.3) 



sup E B I(*)-E^ 2 U 



0<t<l 



1=1 



o P (0 1/2 ) (d^oo) 



Using (B.2), for any 5 G (0, 1), there are a constant c and a positive integer 
Oo such that 

p ( sup J2 B t (*) = su p E B l (*) ) > 1 - ^ 



,0<t<l 



t=l 



if D > 0q • Now 



o< S up 

\t-l/2\<c/^\ i=1 e=1 \*J ) 



sup 

|t-l/2|<c/0 1 /4 



a 8 *(*)-a;{- 



+ o 



t(i — t) - - 

V 1 4 



< y/6 sup 

|t-l/2|<c/0 1 / 4 



a;w-a; 



+ d sup 

0<t<l/2 L 



i(l-i) 



Op 



which follows from the weak convergence in (B.l) and the continuity of 
the limit process applied to the first term on the right-hand side of the 
latter display and obvious reasoning in case of the second term. Relation 
(B.3) is established. To obtain the assertion of Remark 2.1 in the case of 
A(Z>), it suffices therefore to prove asymptotic normality for ^| =1 (1^(1/2) — 
l/4)/(Z)/8) 1 / 2 . This is, however, immediately implied by the CLT. 
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To prove the assertion of Remark 2.1 for 0(c)), we use the relations 
E[Jq Bj{t) dt] = 1/6 and Var(/ Q 1 Bj(t) dt) = 1/45, and the central limit the- 
orem. The proof is now complete. □ 



Proof of Theorem 2.2. It is an immediate consequence of the ergodic 
theorem that, as n — > oo, 



sup 

te[o,9] 



1 \nt\ 

~E^yJ-^[yoy<f] 



sup 

te(0,i] 



\nt\ 



- E yl(yl) T -(*-WyS(yo) T ] 



j=[ne\+l 

By assumption (2.11), it follows that 



op(l), 
o P (l). 



sup 



\nt\ 



; e 



n 



T 

Z j,n z j,n 



1 



LntJ 



^ E 



j=[ne\+l 

Since similar arguments also imply that 



77 

L j=\ne\+l 



l Z j,« Z j,nl 



0P (1). 



sup 



1 



LntJ 



E, T 



j=[n6\+l 



sup 

te(0,i] 



1 



LniJ 



?Z 



E Hn(y*j 



*\T 



j=[ne\+l 



P (1) 



the proof of Theorem 2.2 is complete. □ 



B.2. Properties of the change-point estimator. The purpose of this sub- 
section is to outline a proof that provides us with the relation in (2.13). In 
order to keep the arguments short, we simplify the setting and introduce the 
following technical assumptions. To begin with, we require the pre-change 
sequence (y~j :j<k*) to be as in Assumption 2.2, while the post-change se- 
quence (yj :j>k*) is given by y~j = y| with (yj- :j € Z) as in Assumption 
2.2, thereby omitting the additional noise terms z Jjn here. 

Let Xj = vech[yjyj] and x* = vech[y*(y*) T ] and assume that the thus 
defined sequences satisfy the strong law of large numbers with a rate a n , 
that is, 

■y n 1 n 

— ^(xj — /x) ^ and — E( x i — Z- 4 *) - *0 a,s -> 

where fx = E[x.j] and fx* = E[Xj]. Since we are under Ha, we have that 
A = /i-/x*^0. 
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Theorem B.l. Let k* = [n6] for some 9 G (0, 1). In the above described 
setting, 9 n is strongly consistent for 9 if 

an\ 2 1 n , . 

provided that S n — ► £ mi/i probability one and £ is positive definite. 

If the quantity A n does not depend on the sample size n, then A^S" 1 A > 
converges to a positive constant and the sufficient condition stated in 
Theorem B.l reduces to a n /n — > 0. The latter is essentially satisfied if strong 
laws of large numbers hold for the pre-change and post-change sequences. To 
establish Theorem B.l, we need the following lemma which is given without 
proof. 

Lemma B.l. Let f(t) and g(t) be functions on [0,0] of which f(t) is 
increasing. As long as f(9) — f{6 — 7) > sup t |(/(t)| for some 7 G [0,6], we 
have that 

argmax[/(t) + g(t)] >9-j. 
te[o,0] 

An analogous result can be stated if f(t) and g(t) are functions on [6, 1] of 
which f(t) is decreasing. 

PROOF of Theorem B.l. Observe first that C n = A T E~ 1 A > be- 
cause A 7^ and S" 1 is positive definite by assumption. Next, straight- 
forward computations reveal that 5 n iS~ 1 5 n t = Qi(nt) + Q 2 (nt) + Q 3 (nt), 
where 



Qiint) 



nt 2 (l-9) 2 C n , t€[0,6], 
n6\\-tfL n , te(9,l], 



1 [nt] [nt] n 

Q 2 (nt) = -e(nt) T i)- 1 e(nt), e(nt) = E £ i " — E e J' 

n 3=1 n 3=1 



and 



9) -A T £^e(nt), tG[0,6], 



n 



^Lj-A T ±^s(nt), t€(9,l]. 
Jn 



Q 3 (nt) 



The quadratic form Q\(nt) is clearly increasing in [0,6] and decreasing on 
[6, 1] taking its maximum in 9. It remains therefore to show that Q2{nt) + 
Q- 3 {nt) is small compared to Q\{nt) in the sense of the condition given in 
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Lemma B.l. To this end note that, for vectors a, b 6 M. d and a quadratic 
matrix M G R dxd , \aMb\ < |a||M| B |b|. Thus, for all n > n , 

sup Q2(nt) < — sup ^(rei)^!]" 1 ^ 
te[o,i] n te[o,i] 

< - sup \e{nt)\ 2 \Tr 1 \ E = o(— \ a.s., 
«te[o,i] \nJ 

since < oo by the positive definiteness of S. Further, for all n > uq, 

.„ , 40(1-0). All „ . , /Aa„\ 

sup |Q 3 (nt)| < v ^- lAllS- 1 ^ sup |e(nt)| =o — a.s. 
te[o,i] v n te[o,i] \ V n J 

Since a n / ' yfn — > oo (our sequence satisfies a CLT, therefore the law of large 
numbers cannot hold with a rate smaller than y/n), we have shown that 
sup tg [ ,i] \Q2(nt) + Qi,{nt)\ = o{a 2 l /n) with probability one. Utilizing Lemma B.l 
with fn(t) = Qi(nt) and g n (t) = Q2(nt) + Qz{nt) yields that 

argmax S^E" 1 ^ > 9 - j n , 
te[o,o] 

provided f n (9) -f n {0-ln) = wy n (29-j n )(l-9) 2 C n > sup te [ 0) ^\g n (t)\. A sim- 
ilar argument applies also to the inverval [9, 1]. Choosing 7„ = (a n /n) 2 / C n 
verifies the theorem. □ 



In the special but widely applicable case of weakly dependent sequences 
satisfying the law of the iterated logarithm, one can pick a n = y/n log log re. 
If, in addition, A is also independent of re, then the rate given in (2.13) 
follows. 



B.3. Proofs of Theorems 4.1 4.5. 



Proof of Theorem 4.1. We need to verify the assumptions of The- 
orem A. 2. Note first that the moment assumptions are satisfied. Moreover, 
we can define the approximating random vectors (y^ : j 6 Z) by letting 

Yj 71 ^ = Y^T=o^i £ j-i f° r 3 e ^- Then, condition (2.5) is satisfied, since 
5Zllyo-yo m) || 4 < Ikolk H \C e \<oo 

m>l m>ll>m+l 

by an application of Minkowski's inequality and assumption on the matrices 
{Cg : I > 0). This is the assertion. □ 

Proof of Theorem 4.2. Using the Hadamard product o in the defin- 
ing equations (4.3) and (4.4) allows for proving the assertion componentwise. 
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Let therefore (yj :j£Z) denote a generic coordinate of (y~j : j £ Z). As in the 
statement of Theorem 4.2, we add — if necessary — zero coefficients and as- 
sume that (yj :jeZ) follows the GARCH(r, r) model, r = max{p, g}, 

y j = a j e j , 

r r 

°j = ^ + E a t a ]-e + E PiVj-t- 
1=1 e=i 

The assumption 7c < 1 ensures that (yj :j G Z) has a unique strictly sta- 
tionary solution (see [13, 14]). Berkes, Hormann and Horvath [8] have shown 
that this solution has the representation 

(00 n \ 1/2 

1 + E E nK+^ift-..-^) • 
n=ll<4,.../„<ri=l / 

From this representation and Minkowski's inequality, we obtain that 

(00 n \ 

!+e e niK+^oib 
n=ll<l 1 ,...,e n <ri=l ) 



<oj\\e \\lJ2^c 



n=0 

and consequently E[y^\ < 00 on behalf of 7c < 1. It remains to check con- 
dition (2.5) coordinatewise. In view of (B.4), a natural candidate for the 
approximating variables (y^ : j € Z) is given by the truncated variables 

/ [m/r] „ \l/2 

(B.5) yf > = yfcej 1 + £ £ II K + A,*?-*—*) • 
\ n=l l<ii,...,f„<ri=l / 

Observe that is measurable with respect to the cr-algebra generated by 

the random variables £j,£j_i, . . . ,£j_ m . Using the fact that |yj — yj | 4 < 

|y| — (yj m ' ) ) 2 | 2 , it is now easy to see that (2.5) holds if 7c < 1. The assertion 
of Theorem 4.2 is consequently implied by Theorem A. 2. □ 

Proof of Theorem 4.3. At first, we identify the only candidate for the 
strictly stationary solution iterating the defining volatility equation (4.5). 
Note that, by (4.3), y,- o y„- = aj o crj o e„- o ej, and therefore 

v Q 

o-j o crj = u> + E ° fj-^) + E B e(yj-t ° yj-«) 

r 

(B.6) =u + J2(Ae + B t Ej-i){trj-i o o^) 

4=1 



3G 
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E E I[(M+BeA 



■j-ii- 



Lk=l !<£!,.. .,£ k <ri=l 



(jj. 



where r = max{p, q}, Ej = diag(£j o ej), j € Z. Utilizing the contraction 
condition 7j jQ , < 1, we will show that the expression on the right-hand side 
of (B.6) has a finite || • ||£ jQ ,-norm. To this end, write 



E E iKM+Btity-t!- 

fe=ll<4,...,4<ri=l 



E,a 



oo k 

^e e niK+^oik* 

k=l l<4,...,4<ri=l 



= E E 11^ + B ^\\E,a) 

k=i \e=i J 

oo 

= E < 
k=l 

where the first inequality sign is obtained after an application of the matrix 
norm inequality \MN\e < | A^|^| A^|^, Minkowski's inequality for the L a - 
norm and the fact that (Ej :jgZ) are independent, identically distributed 
random matrices. The proof is complete. □ 



Proof of Theorem 4.4. We need to verify Assumption 2.1. Condi- 
tion (2.3) has been established in Theorem 4.3. It remains to determine 
an approximating sequence (y^ -j £ Z) that satisfies (2.5). To this end, 
introduce 



(m) 

yo 

(m) (m) 



(m) 



E E 



where r = maxjp, q}. Due to the stationarity result given in Theorem 4.3, 
it suffices to prove (2.5) for j = 0. This will be achieved in three steps. We 
first establish a result for the distance between ctq and the approximating 
volatility <x . Observe that, using the definition of the Euclidean norm and 
the L a -norm, we obtain that 



1 0"o 



(m) ,,2 



d 

E 



( a 0,i 



M\2 
0,i J 
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d 

<X]||(o-o,i-^o,i 



i=l 



(B.7) 



^ Ell a o,* - (°"o,i ; 112 
i=i 

^ j II 2 / ("0\2ii 

<d max <r i - (eg j ) I 2 

l<«<a 

j/ || 2 / (m)N2ii2\ 1 / 2 



1/2 



in (m) (™.)n 

= d\\a o a - a y o <x '|| 2 , 

where we have applied that |a — 6| 2 < |a 2 — b 2 \ if a, 6 > to obtain the first 
inequality sign in the second line of the display. We further estimate the 
distance between the Hadamard products of the volatilities ctq o ctq and 



cro™-* o Oq . We have 



I (m) (m) ii 

|<x o cr - <r o <r '\\ 2 



< M 



(B.t 



E E 11(^+^-4--^) 

fc= m +l l<£i,...,4<r i=l 



<h e e nii^+^^oib 

fc=m+l l<4,...,4<ri=l 



= H E En^ + ^oib : 

fc=m+l \£=1 / 

Finally, combining (B.7) with (B.8), we arrive at 

ii (m) ii ii ii (m) ii 

l|yo-y ll 4 < Ikolkllcro -<r ; || 4 



1-7J 



<Vd\\ 



SO 4 



1 -7J 

The right-hand side is summable on account of 7j < 1, thus yielding (2.5) 
and completing the proof of Theorem 4.4. □ 

Proof of Theorem 4.5. It follows immediately from Theorems 4.3 
and 4.4. □ 
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B.4. Proofs of Theorems 4.6 and 4.7. 



Proof of Theorem 4.6. The assertion is an immediate consequence 
of Brandt [15] and Bougerol and Picard [13] . □ 

The proof of Theorem 4.7 requires three auxiliary lemmas and some basic 
results from linear algebra. The first lemma establishes connections between 
the various norms that are in use. 

Lemma B.2. If M is adxd matrix, then it holds that \M\e < V%\ vech[M]|. 

Proof. The inequality follows from the corresponding norm properties. 

□ 

Lemma B.3. Let M and N be two symmetric square matrices, and let 

oo -. oo 

exp(M) = ^— M k and exp(iV) = ^ — N k . 
k=o k - k=o k - 

Then, it holds that: 

(i) exp(M) is positive definite; 

(ii) for every real numbers a and b, exp(aM) exp(bM) = exp((a + b)M); 
(hi) exp(M) exp(— M) = I, where I denotes the identity matrix; 

(iv) (exp(M)) 1 /2 =exp(iM); 

(v) |exp(M + AO-exp(M)| B < \N\ E esxp(\M\ E ) exp(\N\ E ) ■ 

Proof. See, for example, Horn and Johnson [30]. □ 

Lemma B.4. Let X be a positive random variable such that E[exp(tX)] < 
oo if t < r. Then, there is a to £ (0, r) such that E[exp(tX)] < 1 + 2tE[X] 
for all te [0,i ]. 

Proof. The assertion follows from routine arguments. □ 

Proof of Theorem 4.7. To verify Assumption 2.1 note that (2.3) 

follows from Theorem 4.6. It remains to define a suitable sequence (y^ :j G 
Z) that satisfies (2.5). To this end, let 

yo - (-"o ) £ o, 

where the matrix Hq 71 ^ is defined via truncating the strictly stationary so- 
lution at lag m as 

m 

vech[log^ m) -C] = J2 A k F{e- k ^ . . .,e_ fc _,). 

k=0 
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Next, let 



Sq = C + math 



S, 



(m) 



C + math 



^A k F(e- k -i,...,e- k - q ) 

Lfe=0 

m 

5^A*F(e_ fc _i,...,e_ fc _ g ; 

jfc=0 



where math denotes the inverse operator of vech. Then, by assumption on 
C, the matrices Sq and are symmetric. It holds, moreover, that Hq 



exp(So), Hjj m) = exp(5r j ), = exp(±So) and (H^) 



Jm) 



1/2 



r M\l/2 



e*pfrsf>). 

For the last two statements, part (iv) of Lemma B.3 has been applied. Uti- 
lizing the previous statements, we estimate next the matrix norm difference 
between H^ 2 and (Hq^) 1 / 2 as 



\H, 



1/2 



( ff("0)V2| 



|exp(±[S + S£ 



(m) 



< \S - S , ^' m) | £; exp(i|5^ m, | E )exp(i| 1 S"o - S# 



Mi 



:Mi 



1 qM 
2°0 



and iV 



where we have applied part (v) of Lemma B.3 with M 

| (5b — So ) to obtain the inequality sign (additionally dropping a factor \ 

on the right-hand side) . Since is independent of Sq and Sq"^ , we conclude 
further that 



|yo 



Mm 

yo U 

< I eo \H { 



1/2 



r rrM\l/2|| 

l-"o > Wea 



(B.< 



<||6 ||4|||So-Sr ) | i3 exp(I|4 m} |i ? )exp(I|5o-4 
< lleolkllSb - 5 , ^ m) || £ ; )8/3 ||exp(i|S'^ m) | £; )|| 4Q , 



?Mi 



x ||exp(||Sb-^o |b)II8/3- 

The latter inequality sign follows from an application of Holder's inequal- 
ity with a,j3 > 1 such that 1/a + 1/(5 = 1. In the following, we investigate 
the norms in (B.9). By assumption, ||eo||4 < oo. Lemma B.2, Minkowski's in- 
equality and the existence of a moment generating function for F(s-i, . . . , e_ g 
imply that 



||5 -5, 



) V^<v^||vech[S -Sf ) ]|| 8/3 



< V2 



£ \A\ k E \F(e_ k „ u 

k=m+l 



Hi3 
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<V2 £ |A|*||F( e _ fc _ 1) ... ) e_ fc _,)|| 8i9 

fc=m+l 

<C|A|£<oo, 

where C > is a constant, on account of \A\e < 1 by assumption. By similar 
arguments, we continue estimating 



expQ|So-Sj" } |s) 

exp^-t=|vech[S - 5"J m) ]|^ 



< 



< 



exp 



< 



y/2 
1 



fc=m+l 



8/3 



exp H= X! |^|ll^(e-fe-i,---,e-fe- 9 ) 
X(m,/3). 



8d 



Observe that (F(e_k-i, ■ ■ ■ > £-k-q) '■ k > 0) is a sequence of (/-dependent ran- 
dom variables. Therefore, each of the sequences . . . , e_k-q) '■ k = 
z^modg), z/ = 0, 1, . . . ,q — 1, consists of independent random variables. Defin- 
ing 

I u {m)= \A\ k E \F(£_ k _ u ...,e- k „ q )\, 

m<k=u mod q 

we obtain that 



v = 0,l,...,q- 1, 



J(m,/3) 



J5 



TT f 8/3 r/ ^ 
11 exp^-^=^(m) 



,i/=0 



< 

\i/=0 



6XP V \72 



V(8/3) 



1/(8/3?) 



It remains to estimate the expectations in the latter equation. We do so first 
assuming that m is sufficiently large. In this case, the fact that /^(m) is a 
sum of independent random variables and Lemma B.4 imply that there is a 
constant c such that 



E 



exp 



(l 4(m) 



n e 

m<k=v mod q 



exp 



On*.-.,...... 
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< [] (l + c|A|iE[|F(e_i,...,e- 

m<k=u mod q 

< [] exp(c\A\ k E E[\F(e^,...,e. 

m<k=v mod q 

<exp(c|A|£), 



if m > itlq for a suitable constant c > 0. If m < mo, we set .Hq 177 ^ = 0. Then 
the proof is complete if we show that 



sup 

Kra<oo 



'( \\ S m) \E 



< oo, 



la 



where S^ ^ = Sq. This, however, follows from similar arguments as before as 



long as the moment generating function E[ex.p(t\F(s-i, . . . ,S- q )\)] 
for t > Aaq/y/2. Now pick a close to 1 to obtain the assertion. □ 



is finite 



REFERENCES 

[1] Andreou, E. and Ghysels, E. (2002). Detecting multiple breaks in financial market 
volatility dynamics. J. Appl. Econometrics 17 579-600. 

[2] Andreou, E. and Ghysels, E. (2003). Tests for breaks in the conditional co- 
movements of asset returns. Statist. Sinica 13 1045-1073. MR2026061 

[3] Andrews, D. W. K. (1984). Nonstrong mixing autoregressive processes. J. Appl. 
Probab. 21 930-934. MR0766830 

[4] Andrews, D. W. K. (1991). Heteroskedasticity and autocorrelation-consistent co- 
variance matrix estimation. Econometrica 59 817-858. MR1 106513 

[5] Aue, A., Berkes, I. and Horvath, L. (2006). Strong approximation for the sums 
of squares of augmented GARCH sequences. Bernoulli 12 583-608. MR2248229 

[6] Bai, J. and Perron, P. (1998). Estimating and testing linear models with multiple 
structural changes. Econometrica 66 47-78. MR1616121 

[7] Bauwens, L., Laurent, S. and Rombouts, J. V. K. (2006). Multivariate GARCH 
models: A survey. J. Appl. Econometrics 21 79-109. MR2225523 

[8] Berkes, I., Hormann, S. and Horvath, L. (2008). The functional central limit 
theorem for a family of GARCH observations with applications. Statist. Probab. 
Lett. 78 2725-2730. MR2465114 

[9] Berkes, I., Horvath, L. and Kokoszka, P. (2003). GARCH processes: Structure 
and estimation. Bernoulli 9 201-227. MR1997027 
[10] Billingsley, P. (1968). Convergence of Probability Measures. Wiley, New York. 
MR0233396 

[11] Bollerslev, T. (1986). Generalized autoregressive conditional heteroskedasticity. 
J. Econometrics 31 307-327. MR0853051 

[12] Bollerslev, T. (1990). Modeling the coherence in short-run nominal exchange rates: 
A multivariate generalized ARCH model. Rev. Econom. Statist. 74 498-505. 

[13] Bougerol, P. and Picard, N. (1992). Strict stationarity of generalized autoregres- 
sive processes. Ann. Probab. 20 1714-1730. MR1188039 

[14] Bougerol, P. and Picard, N. (1992). Stationarity of GARCH processes and of 
some nonnegative time series. J. Econometrics 52 115-127. MR1 165646 



42 



AUE, HORMANN, HORVATH AND REIMHERR 



[15] Brandt, A. (1986). The stochastic equation Y n+1 = A n Y„ + B n with stationary 
coefficients. Adv. in Appl. Probab. 18 211-220. MR0827336 

[16] Brockwell, P. J. and Davis, R. A. (1991). Time Series: Theory and Methods, 2nd 
ed. Springer, New York. MR1093459 

[17] CARRASCO, M. and Chen, X. (2002). Mixing and moment properties of vari- 
ous GARCH and stochastic volatility models. Econometric Theory 18 17-39. 
MR1885348 

[18] CSORGO, M. and Horvath, L. (1997). Limit Theorems in Change-Point Analysis. 
Wiley, Chichester. 

[19] Davidson, J. (1994). Stochastic Limit Theory. Oxford Univ. Press, Oxford. 
MR1430804 

[20] Davis, R. A., Lee, T. C. M. and Rodriguez- Yam, G. (2008). Break detection for a 
class of nonlinear time series models. J. Time Ser. Anal. 29 834-867. MR2450899 

[21] Davis, R. A. and MlKOSCH, T. (1998). The sample autocorrelations of heavy-tailed 
processes with applications to ARCH. Ann. Statist. 26 2049-2080. MR1673289 

[22] Doukhan, P. and Louhichi, S. (1999). A new weak dependence condition and appli- 
cation to moment inequalities. Stochastic Process. Appl. 84 313-343. MR1719345 

[23] Engle, R. F. (1982). Autoregressive conditional heteroskedasticity with estimates 
of the variance of U.K. inflation. Econometrica 50 987-1008. MR0666121 

[24] Engle, R. F., Ng, V. K. and Rothschild, M. (1990). Asset pricing with a factor 
ARCH covariance structure: Empirical estimates for treasury bills. J. Econo- 
metrics 45 213-238. 

[25] Erdelyi, A. (1954). Tables of Integral Transforms 1. McGraw-Hill, New York. 

[26] Gombay, E., Horvath, L. and Huskova, M. (1996). Estimators and tests for 
change in the variance. Statistics and Decisions 14 145-159. MR1406521 

[27] Hall, P. and Yao, Q. (2003). Inference in ARCH and GARCH models with heavy- 
tailed errors. Econometrica 71 285-317. MR1956860 

[28] He, C. and Terasvirta, T. (1999). Fourth moment structure of the Garch(p, q) 
process. Econometric Theory 15 824-846. 

[29] Hormann, S. (2008). Augmented GARCH sequences: Dependence structure and 
asymptotics. Bernoulli 14 543-561. 

[30] Horn, R. A. and Johnson, C. R. (1991). Topics in Matrix Analysis. Cambridge 
Univ. Press, Cambridge. MR1091716 

[31] Hsu, D. A. (1979). Detecting shifts of parameter in gamma sequences with appli- 
cations to stock price and air traffic flow analysis. J. Amer. Statist. Assoc. 74 
31-40. 

[32] Ibragimov, I. A. (1962). Some limit theorems for stationary processes. Theory 
Probab. Appl. 7 349-382. MR0148125 

[33] Inclan, C. and Tiao, G. C. (1994). Use of cumulative sums of squares for retro- 
spective detection of changes of variance. J. Amer. Statist. Assoc. 89 913-923. 
MR1294735 

[34] Jeantheau, T. (1998). Strong consistency of estimators for multivariate ARCH mod- 
els. Econometric Theory 14 70-86. MR1613694 

[35] Kawakatsu, H. (2006). Matrix exponential GARCH. J. Econometrics 134 95-128. 
MR2328317 

[36] Kiefer, J. (1959). if-sample analogues of the Kolmogorov-Smirnov and Cramer-v. 

Mises tests. Ann. Math. Statist. 30 420-447. MR0102882 
[37] Kokoszka, P. and Leipus, R. (2000). Change-point estimation in ARCH models. 

Bernoulli 6 513-539. MR1762558 



MULTIVARIATE COVARIANCE BREAK DETECTION 



4.3 



[38] Lanne, M. and Saikkonen, P. (2007). A multivariate generalized orthogonal factor 
GARCH model. J. Bus. Econom. Statist. 25 61-75. MR2338871 

[39] Li, W. K. (2004). Diagnostic Checks in Time Series. Chapman & Hall, Boca Raton, 
FL. 

[40] MlKOSCH, T. (2003). Modeling dependence and tails of financial time series. In: Ex- 
treme Values in Finance, Telecommunications, and the Environment (B. Finken- 
staedt and H. Rootzen, eds.) 185-286. Chapman & Hall, Boca Raton, FL. 

[41] MlKOSCH, T. and Straumann, D. (2006). Stable limits of martingale transforms 
with application to the estimation of GARCH parameters. Ann. Statist. 34 493- 
522. MR2275251 

[42] Nelson, D. B. (1991). Conditional heteroskedasticity in asset returns: A new ap- 
proach. Econometrica 59 347-370. MR1097532 

[43] Nze, P. A. and Doukhan, P. (2004). Weak dependence: Models and applications 
to econometrics. Econometric Theory 20 995-1045. MR2101950 

[44] Petrov, V. V. (1995). Limit Theorems of Probability Theory. Oxford Univ. Press, 
Oxford. MR1353441 

[45] Potscher, B. M. and Prucha, I. R. (1997). Dynamic Nonlinear Econometric Mod- 
els, Asymptotic Theory. Springer, New York. MR1468737 

[46] Serban, M., Brockwell, A., Lehoczky, J. and Srivastava, S. (2007). Modelling 
the dynamic dependence structure in multivariate financial time series. J. Time 
Ser. Anal. 28 763-782. MR2395913 

[47] Silvennoinen, A. and Terasvirta, T. (2008). Multivariate GARCH models. In: 
Handbook of Financial Time Series (T. G. Andersen, R. A. Davis, J. -P. Kreiss 
and T. Mikosch, eds.) 201-229. Springer, New York. 

[48] Tsay, R. S. (2005). Analysis of Financial Time Series, 2nd ed. Wiley, New York. 
MR2162112 

[49] Vrontos, I. D., Dellaportas, P. and Politis, D. N. (2003). A full- factor multi- 
variate GARCH model. Econom. J. 6 312-334. MR2028238 

[50] Wu, W. B. (2005). Nonlinear system theory: Another look at dependence. Proc. 
Natl. Acad. Sci. USA 102 14150-14154. MR2172215 



A. Aue 

Department of Statistics 
University of California 
One Shields Avenue 
Davis, California 95616 
USA 

E-MAIL: alexauc@wald.ucdavis.edu 



s. hormann 
l. horvath 

Department of Mathematics 

University of Utah 

155 South 1440 East 

Salt Lake City, Utah 84112-0090 

USA 

E-MAIL: hormann@math.utah.edu 
horvath@math.utah.edu 



M. Reimherr 

Department of Statistics 

University of Chicago 

5734 Street University Avenue 

Chicago, Illinois 60637 

USA 

E-MAIL: mreimherr@uchicago.edu 



